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ABSTRACT 

We present results from integral field spectroscopy with the Potsdam Multi- 
Aperture Spectrograph at the 3.5m telescope at Calar Alto Observatory of the intense 
star-forming region [HL90] 111 at the center of the starburst galaxy IC 10. We have 
obtained maps with a spatial sampling of 1" x 1" = 3.9 pc x 3.9 pc of different emis- 
sion lines and analyzed the extinction, physical conditions, nature of the ionization, 
and chemical abundances of the ionized gas, as well determined locally the age of the 
most recent star-formation event. By defining several apertures, we study the main 
integrated properties of some regions within [HL90] 111. Two contiguous spaxels show 
an unambiguous detection of the broad He II A4686 emission line, this feature seems 
to be produced by a single WNL star. We also report a probable N and He enrichment 
in the precise spaxels where the WR features are detected. The enrichment pattern is 
roughly consistent with that expected for the pollution of the ejecta of a single or a 
very small number of WR stars. Furthermore, this chemical pollution is very localized 
(~ 2" ~ 7.8 pc) and it should be difficult to detect in star-forming galaxies beyond 
the Local Volume. We also discuss the use of the most-common empirical calibrations 
to estimate the oxygen abundances of the ionized gas in nearby galaxies from 2D 
spectroscopic data. The ionization degree of the gas plays an important role when 
applying these empirical methods, as they tend to give lower oxygen abundances with 
increasing ionization degree. 
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1 INTRODUCTION 

Dwarf and irregular galaxies are ideal laboratories to study 
the interaction between massive stars and the interstellar 
medium (ISM). Indeed, massive stars do not only provide 
the energetic radiation that ionizes the H I clouds into H II 
regions, but they also inject mechanical energy into the sur- 
rounding medium, both via strong stellar winds and super- 
nova explosions, which create shells and bubbles and may 
induce the birth of new stars. The effects of the massive 
stars in the ISM are especially strong when they are gath- 
ered together in large clusters; their actions can even af- 



fect the nearby inter-galactic medium (IGM). Strong star- 
forming regions, which may be found in galaxies covering 
a wide range of masses, may enrich both the ISM and the 
IGM. Actually, it is still unclear whether the chemical en- 
richment provided by massive stars is an instantaneous event 
or if the new synthesized heavy ele ments must cool be- 
fore becoming part of the ISM (e . g., Tenorio-TaEle 19961; 
Kobulnickv fc Skilmianlll996l, Il997l, 1 19981; iKobulnickv et al 



1997 



Henry. Edmunds fc Koppenll200(]| ; Ivan Zee fc Havne; 



20061 ; lLopez-Sanchez fc Estebanll2010bl 'l. All these processes 
have a large effect on the formation and evolution of galax- 
ies. 



* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by 
the Max-Plank Institut fur Astronomie and the Instituto de As- 
trofisica de Andalucfa (CSIC). 

f Visiting Astronomer at the Instituto de Astrofisica de Canarias. 
E-mail: Angel. Lopez-Sanchez@csiro.au 



Nearby dwarf galaxies have the advantage that they 
can be observed at high spatial resolution and even, in 
some cases, their stars are resolved. This allows a de- 
tailed analysis of the interplay between the neutral and 
ionized components of the ISM and the massive stars 
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Figure 1. {Left) False-color image of IC 10 obtained combining data provided bv lGil de Paz, Madore fc Pevunoval i2003h in B (blue), R 
(green) and Ha (red) filters. The FoV is 6.8' X 7.4', which covers the majority of the BCDG as it is seen in optical light. Notice the long 
and extensive ionized gas emission throughout all the galaxy. The main star-forming region is shown by a white box. (Right) Enlargement 
of th e main star-forming region within IC 10 showing the pure Ha emission (continuum-subtracted) of the galaxy ( Hunter & El megreer] 
2004). Its FoV is 57"x57". The blue box correspon ds to the field of vie w of the IFU used with PMAS, which has a size of 16"xl6" and 
covers the intense star-forming region [HL90] 111c l|Hodge fc Leelll990l . see their Fig. 2b). The ident ification of the different H II regions 
follows the names given by these authors. The positions of the individual WR stars known in IC 10 llCrowther et al]|2003t) are shown by 
stars. Three of them, [MAC92] 24 A, B and C (yellow stars), are located at the center of our observed area. 



(e.g.. iThurow fc Wilcotj 2005: van Evmeren et al ] 12009a b. 



20101; iMonreal-Ibero et alj |2010| )~ and constrains the kine- 
matical and chemical evolution of the galaxies. Recently, 
our group confirmed a localized chemical enrichment in 
the nearby (4 Mpc) blue compact dw arf galaxy (BCDG) 
NGC 5253 (lLopez-Sanchez et al.ll2007l ). The nitrogen pol- 
lution (and probably also helium) found in a particu- 
lar area of the galaxy seems to be a consequence of 
the ejecta of Wolf-Rayet (WR) stars. Evidences of pos- 
sible nitrogen pollution has been only reported in some 
few starburst s, but the majority of them also showing 
WR features riKobulnickv et al.lll997l; iPustilnik et alj|2004 
Ijames et al.l 20091 ; lLopez-Sanchez fc Estebanl l2010al lbf). 
However, NGC 5253 remains as the only known star- 
forming galaxy where a localized che mical enrich- 
ment has been defini tively established (I Walsh fc 



19891; iKobulnickv et al" 



Monreal-Ibero et al.ll2010 ) 



19971 ; lLopez-Sanchez et al.1 1 200 



The advent of the 2D spectroscopy allows to ex- 
plore this issue more efficiently, as this technique permits 
to record simultaneously the optical/near-infrared (NIR) 
spectra of regions within the galaxies without the limita- 
tions of long-slit observations. In particular, spectroscopy 
of dwarf, star-forming galaxies using integral field units 
(IFU) provides a very powerful tool to constrain the phys- 
ical, chemical and kinematic al properties of the ionized 
gas within these objects (e.g.. [Garcia- Lorenzo^ et aT1l2008l; 
iLagos et alj|2009l ; ICairos et all l2009al lbl l201fj ; Ijames et al.1 



2009, |2010D . to look for the places where massive stars 
are located, to study the underlying stellar population 
and its kinematics via the analysis of the continuum and 
the absorption lines, and finally to combine all the re- 
sults to understan d the feedback between massive stars 
and th e ISM (e.g.. IThurow fc Wilcotj 120051; iKehrig et al.1 



| 200Sl; lAlonso-Herrero et alj 120091 ; iBordalo et all 1200a 
Monreal-Ibero et al.ll2010l ). Here, we present an optical 2D 



spectroscopical analysis of a particularly interesting and in- 
tense star-forming region within the nearest starburst, the 
galaxy IC 10. Although our main aim is to look for a local- 
ized chemical enrichment within the H 11 region as a conse- 
quence of the massive stars, we will also analyze the extinc- 
tion, ionization structure, chemical abundances and age of 
the most recent star-formation burst in this particular area 
of IC 10. 

Sometimes defined as an irregular galaxy, IC 10 
(Fig. [1} actually is the only known st arburst with in the 
Local Group. Because of it s low mass dMateol 19981). high 



star-f ormation activity (i.e.. lHodge fc Leej|l9 90; Lero v et al 



l200fj) l ow metall ic ity -1 2 +log(Q/H)=8.26. ISkillman et al. 
jimt: iGarnettJ l|l990l ); iMagrini fc Goncalvej (|2009l )-. 



its very extended, clumpy, and pec uliar distribution of 
neutr a l gas (Juette et al., in prep; [Shostak fc Skillmanl 
1 19891 ; IWilcots fc Miller! [ 1998), high far-infrared emis- 
sion ( Melisse fc Israel Il994l), non-thermal radio- 



continuum emission ( Yang fc Skillmanl [l993), and the 
relative importance of its young stellar population 
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(i.e., ISakai. Madore fe Freedman 1999|; Borissova et ah! 
l200d ; iHunterl l200ll ; IVacca. Sheehv fc Grahaml 120071 ) 
which includes a large numb er of WR stars 
(IMassev. Armandroff fc Conti 1 19921 ; iMassev fc Johnson! 
1 19981 ; IMassev fc Holmes! 12003). IC 10 is usually classified 
as a BCDG (|Richer et al] l200ll ). The high carbon-rich 
WR stars to nitrogen-rich WR stars ratio in IC 10 (the 
WC/ WN ratio), which is u nusually high for its metal- 
licity (|Crowther et al.l [2003), suggests that the galaxy 
has experience a brief but intense galaxy-wide burst of 
star-formation within the last 10 Myr. The starburst was 
probably trigg ered by gas falling in from the large, exte nded 
H i reservoir jsaito et alJll993 ; IWilcots fc Millerlll99Sl ). 

The determination of the distance to IC 10 has been 
challenging because of its proxi mity to the Galactic plane 
(I = 1 19.0° b = -3.3°) Recen tly, IVacca, Sheehv fc Grahaml 
(|2007l ) and lSanna et alJ (|200ST ) used the tip of the red giant 
branch (TRGB) to compute a distance of 0.79 and 0.81 Mpc, 
respectively, to this galaxy. Throughout this paper, we will 
consider that IC 10 is located at 0.80 Mpc, and hence the 
spatial scale is 1" = 3.88 pc. 

The structure of this paper is the following. In we 
describe the observations and the reduction procedure. ij3] 
presents the measurements of the emission lines and how the 
reddening correction was performed, jjjeontains the analysis 
of the ionized gas (extinction, physical conditions, nature 
of the ionization and chemical abundances) as well as the 
localization of the WR stars and the determination of the 
age of the most recent star- formation event. ij5] discusses the 
two main results found in our analysis: the detection of a 
probable N and He enrichment in the precise place where 
the WR features are detected and a warning about the use 
of empirical calibrations to estimate the oxygen abundances 
of the ionized gas in nearby galaxies from 2D spectroscopic 
data. We summarize our conclusions in 36] 



2 OBSERVATIONS AND DATA REDUCTION 



Ha emission ma 



The massive star-forming region [HL 901 111 dHodge fc Led 
Il990l ) , also known as [LMR79] #1 jLequeux et al] 1 19791 ). 
within the starburst galaxy IC 10 was observed on 2007 
October 14 at Calar Alto Observatory (Almerfa, Spain), us- 
ing the 3.5m Telescop e with the Potsdam Multi-Apertur e 
Spectrometer (PMAS. [Roth et all I200H iKelz et~alTl2006t ). 
PMAS is an integral field spectrograph, with a lens array of 
16x16 square elements, connected to a bundle of fiber op- 
tics, whose 256 fibers are re-arranged to form a pseudoslit 
in the focal plane of the spectrograph. We used the stan- 
dard lens array IFU of 16x16 arcsec 2 field of view (FoV) 
with a sampling of 1" which, at the distant to IC 10, pro- 
vides a spatial sampling of 3.88 pc/arcsec and a FoV of 
62.1x62.1 pc 2 . Most of the optical range was covered with 
the V600 grating using two grating rotator angles: —12°, 
covering from 3500 to 5070 A; and —68°, covering from 
5700 to 7150 A. The effective spectral resolution was 3.6 
A. The blue and red spectra have a total integration time 
of 1800 seconds in both cases. The position of the IFU cov- 
ering the massive star- forming region [HL90] 111 is showed 
in Figure [T] Its center has the coordinates a=00 h 20 m 2719, 
5=59° 17' 35.9". The reconstruction of the Ha map using 
our 2D PMAS spectroscopy (Fig. [2]) very nicely matches 




a R.A. [ arcsec ] 



Figure 2. Ha flux map of the intense star-forming region 
[HL90] 111 l|Hodge fc Ledll99(t ) from our PMAS data, which has 
a size of 16x16 pixels (each pixel has a size of 1" X 1"). The flux 
is in units of erg cm -2 s — 1 and it is not corrected for extinction. 
The positions of the broad He II A4686 emission feature, which 
indicates the presence of WR stars, are shown by a star (strong 
blue WR bump) and a cross (detection of the blue WR bump). 
The o symbol indicates the position of a tentative detection of 
the nebular He II A4686 emission line. The four different regions 
used to perform a more detailed analysis of the ionized gas are 
also shown. 



the Ha image provided bv lHunter fc Elmegreenl (|2004l ) and 
showed in the right panel of Fig[T] The PMAS data showed 
two spaxels which poor light transmission, [1,15] and [8,16], 
which we will not consider in our analysis. 

The giant H 11 region [HL90] 111c (see Fig. [TJ harbors 
one of the brigh test WR stars in the galaxy, that was named 
[MAC92] 24 by IMassev. Armandroff fc Contil (I1992T) . Origi- 
nally identified as a WN star bv lRicher et al.l l|200ll ) because 
of the detection of the broad He 11 A4686 emission and the 



this object is a blend of several stars (IMassev & Holmes 


2002 


; ICrowther et al. 


2003|; Vacca, Sheehv & Graham 


2007 


). In particular, 


Crowther et al.l (120031) noted that 



there were three closely spaced sources, [MAC92] 24A, B, 
and C (see Fig[T]), al l located within 1" — 2" at the so uth of 
[HL90] 111c. Later, IVacca. Sheehv fc Grahaml (|2007l ) used 
ground-based laser guide star adaptive optics at the Keck II 
telescope to resolve [MAC92] 24A in at least three bright, 
blue sources and [MAC92] 24B into two blue sources, sug- 
gesting that there are at least four robust WN candidates in 
the field. This hypothesis should be confirmed using optical 
spectroscopy. And, indeed, our 2D PMAS spectra indicate 
the detection of the blue WR bump precisely at the two 
spaxels where the WN star candidates [MAC92] 24 B,C 
should be located. Their position is indicated with a star 
(spaxel [11,8], which shows a strong blue WR bump) and a 
cross (spaxel [11,9], where the WR bump is also detected 
but not as intense as in the previous position) in Fig [21 as 
we discuss below. The maximum of Ha emission provided 
by our 2D spectroscopic data is located just 1" — 2" at the 
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north of the position of the WR stars, and it coincides with 
th e center of the intense r adio-continuum source reported 
bv lYang fc SkillmarJ (|l993h . [YS93] 001744.6+590101. 

Calibration images were obtained during the night: 
arc lamps for the wavelength calibration and a contin- 
uum lamp needed to extract the 256 individual spectra on 
the CCD. Observations of the spectrophotometri c standard 
stars BD +28°4211, Feige 110 and G 191-B2B i|Okelll990l ) 
were used for flux calibration. The absolute error of this cal- 
ibration was of the order of 5%. Night was photometric and 
the typical seeing during the observations was 1". 

The data were reduced using the iRAlj]] reduction pack- 
age SPECRED. After bias subtraction, spectra were traced on 
the continuum lamp exposure obtained before each science 
exposure, and wavelength calibrated using HgNe arc lamp. 
The continuum lamp and sky flats were used to determine 
the response of the instrument for each fiber and wavelength. 
Finally, for the standard stars we have co-added the spec- 
tra of the central fibers and compared with the tabulated 
one-dimensional spectra. 

The effect of the differential atmospheric refraction 
(DAR) is negligible because of both the low airmass (be- 
tween 1.08 and 1.11) in which IC 10 was observed and the 
short exposure time (~1 hour). 
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Figure 3. Example of the optical spectrum obtained at the cen- 
ter of the BCDG IC 10 using a single fiber of our PMAS data. 
In particular, we plot a spectrum with a high S/N, which corre- 
sponds to the maximum of Ha emission (spaxel [11,12]). The most 
important emission lines are labeled. The spectrum has been cor- 
rected for radial velocity, but not for extinction or sky emission. 
In particular, note the Hg I contamination at A4359. 



3 LINE MEASUREMENTS 

We measured all the emission lines detected in the spectra: 
hydrogen Balmer lines (Ha, H/3, H/y, H<5, and, in some cases, 
H8 and H9), helium lines (He I A4471, A5875, A6678 and 
A7065) and the typical collisional excited lines of metallic 
elements: [O n] AA3724,29 (blended), [O in] AA4959,5007, 
[N II] AA6548,83, [S II] AA6716,31, [Ar III] A7136 and 
[Nc in] A3869 and A3967 (blended with H7). Figure [3] shows 
an example of the spectra provided by our PMAS observa- 
tions in which all important emission lines are labelled. It 
corresponds to the spectrum obtained in the region with the 
maximum of Ha emission (spaxel [11,12]). 

Line fluxes were measured integrating all the emission 
between the limits of the line and over a local adjacent con- 
tinuum. All these measurements were made with the SPLOT 
routine of iraf and using our o wn scripts to achieve an au - 
tomatic process (see details in iMesa-Delgado et al.l F2009). 
However, because of the faintness of many of the emission 
lines, a visual inspection of the spectra was needed to get 
a proper estimation of the adjacent continuum and the line 
flux in these cases. We want to emphasize here the impor- 
tance of this manual step because, although automatic pro- 
cedures usually work well when measuring high signal-to- 
noise emission lines, they could give a wrong estimate of the 
flux of faint emission lines. 

The errors associated with the line flux measurements 
were determined following the error cri t erion described in 
IMesa-Delgado. Esteban fc Garcfa-Roiasl i|2008l ). The final 
error of a line was the quadratic sum of the error in its 
flux measurement and the error in flux calibration. In order 
to avoid spurious weak line measurements, we imposed three 



1 IRAF is distributed by NOAO which is operated by AURA Inc.. 
under cooperative agreement with NSF. 



criteria to discriminate between real features and noise: 1) 
line intensity peak over 2.5 times the sigma of the contin- 
uum; 2) FWHM(A) > 1.5xFWHM(H i) and FWHM(A) < 
FWHM(H i)/1.5; and 3) F(A) > 0.001xF(H/3). 

All line fluxes of a given spectra have been cor- 
rected for extinction using the theoretic al Ha/H/3 ratio ex- 
pected for case B recombination given bv lStorev fc Hummer! 
( 1995) assuming the physical conditions of T e = 10000 K and 
n e = 100 cm -3 , the expected values for region [HL90] 111 
in IC 10 (see JO), [Ha/H/3] thc o=2.86. We used the 
reddening function, /(A), norma lized to H/3 derived by 
ICardelli. Clayton fc Mathil (|l989h to correct for reddening 
each emission line, following the expression 



7(A) _ F(X) 



/(H/3) F(H/3) 



x 10 



c(H/3)x/(A) 



(1) 



being /(A) and F(X) the real and the observed fluxes, respec- 
tively, of the emission line and c(H/3) the reddening coeffi- 
cient derived from the Ha/H/3 ratio. As we describe below 
(and as it should be expected for IC 10 because of its very 
low Galactic latitude) the reddening correction is fundamen- 
tal to get a proper estimation of the line fluxes, and indeed 
the derived c(H/3) almost always provides values higher than 
1.0 dex. 

Unfortunately, the Hg I A4358 telluric emission was 
completely blended with the faint [O ill] A4363 auroral line, 
which has a redshifted wavelength of A4359 in IC 10. Using 
the regions with the faintest ionized gas emission (AP 2, see 
below), we have estimated an average value of the Hg I A4358 
line flux of 1.8x 10 -16 erg cm -2 s _1 for each spaxel. In these 
resolution elements, and assuming an electron temperature 
of T e ~10000 K, we estimate that the [O ill] A4363 line flux is 
~80-100 times fainter than that measured for the Hg I A4358 
line flux and hence completely negligible. However, the con- 
tribution of the [O ill] A4363 line to the Hg I A4358 line 
flux is important in the brightest regions of the nebula. In 
these cases, and applying the same method, we estimate that 
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Table 1. Dereddened line intensity ratios with respect to J(H/3)=100 for the integrated spectrum of the intense star-forming region 
[HL90] 111 at the starburst galaxy IC 10 and the regions we specifically analyse within it. Other interesting properties derived from 
these data are also listed. Colons indicate measurements with an uncertainty larger than 40%. Double colons indicate a very tentative 
measurement. Three dots indicate that the emission line has not been detected. 
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88.3+3.6 


10.0+1.3 


57.7+3.6 


11.3+1.1 


14.9+1.3 


167+8 


-W H7 [A] 






13.81+0.12 


28.8+2.6 


3.3+0.9 


19.3+1.9 


3.2+0.6 


5.7+2.3 


38.4+1.1 


-W HS [A] 






6.8+1.1 


11.4+1.3 




9.9+1.2 




3.02: 


17.7+4.5 


A Ka [mag] 




2.49+0.09 


2.51+0.14 


2.51+0.16 


2.35+0.07 


2.18+0.07 


2.02+0.07 


2.61+0.04 


/(Ha) [10- 


-14 —1 
erg s cm 




1547+45 


485+19 


14.5+0.9 


70.0+2.4 


22.3+1.1 


7.55+0.42 


36.1+2.4 


L(Ha) [10 36 ergs" 1 ] 




1185+64 


371+34 


11.1+1.3 


53.6+4.5 


17.1+1.5 


5.78+0.54 


27.6+2.4 


Mhii [M ] 




17600+950 


5510+510 


165+19 


796+67 


254+22 


86+8 


410+36 


M Hn /area [Mq pc~ 2 ] 




4.60+0.25 


15.3+1.4 


1.10+0.13 


4.81+0.40 


8.43+0.73 


5.71+0.53 


27.2+2.4 



F([0 in] A4363)~0.4-0.6xF(Hg I A4358). The contribution 
of the auroral line is even more important when integrating 
the line fluxes of several fibers, as we do in several defined 
apertures. Hence, we will use the derived average value of 
the Hg I A4358 line flux to get a estimate of the flux of the 
[O ill] A4363 line within these areas. 



4 RESULTS 

The reconstruction of the Ha map using our 2D PMAS spec- 
troscopy is shown in Fig. [2] This map cor r espon ds to the 
area named as [HL90] 111 bv lHodge fc Lee] (|l990h . We dis- 
tinguish a very strong star-forming region at the NW (region 
c), which has a size of ~ 10" x 6". Other fainter H II regions 
are located at the E (knot d) and in the SW corner (knot a), 
this last one is only partially covered by our FoV (see Fig.[T|. 
An important feature is a hole located at the south of the 
brightest star-forming region. Indeed, there is an arc-like fil- 



ament of ionized gas (knot e following iHodge fc Ledll99ut ) 
connecting the SW tip of the brightest star- forming area (c) 
with the H II region located at the east (d). The minimum 
of Ha emission (the hole) is located between these three 
structures. 

We have integrated the spectra within several regions 
to analyse in detail the properties of the ionized gas and 
to maximize the signal-to-noise (S/N) ratio of the ID- 
spectra. Figure [2] shows the position of these regions: AP 1 
(which corresponds to the brightest star-forming region, 
[HL90] 111c), AP 2 (which corresponds to the hole of Ha 
emission), AP 3 (which corresponds to [HL90] Hid) and 
AP 4 (which integrates the spectra of the two spaxels where 
the WR features are detected). For these four regions we 
measured all the line fluxes following the same procedure 
explained before. We also carefully analysed two individ- 
ual fibers: #90 (spaxel [11,10]), where a tentative detection 
of the nebular He II A4686 line is found, and #92 (spaxel 
[11,13]), which corresponds to the maximum in the Ha emis- 
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c(Hjg) (Extinction mop) _+ Hq contours 




Density mop from [S II] X6716 / [S II] X6731 [cm" 3 




Figure 4. {Left) Map of the reddenning coefficient, c(H/3), derived from the Hct/H/3 ratio in the field observed in IC 10. (Right) Density 
map derived from the [S II] A6716/A6731 ratio. The Ha contours are overlaid in both maps. Symbols are the same that in Fig. [2] 



sion and hence to the fiber where emission lines have the best 
S/N ratio. 

Because of the stacking of the spectra, some emis- 
sion lines that were not easily observed in the individ- 
ual spaxels (i.e., H9, H8, H<5, [Ne ill] A3869, [Ne III] 
A3967+He, He I A4471) are now clearly distinguished. That 
allowed us to perform a more accurate determination of 
the reddening coefficient within these areas. We consid- 
ered all available H I Balmer line ratios (Ha/H/3, Rj/R/3, 
R8/R/3, and also H8/H/3 in the brightest region) and ap- 
plied the formalism explai n ed in lL6pez-Sanche3 <|2006n and 
lLopez-Sanchez Sz Estebanl |2009l ) to simultaneously derive 
c(H/3) and the equivalent width of the absorption in the H I 
lines, Wabs, which we assumed to be same for all H I Balmer 
lines. Table [1] shows the dereddened line intensity ratios and 
their associated errors, as well as the adopted /(A) of each 
emission line, measured for the total integrated spectrum, 
the four regions we have carefully analyzed, and fibers #90 
and #92. In this table we also include other important quan- 
tities such as the size of the extracted aperture, the observed 
H/3 flux (uncorrected for extinction), the adopted values of 
c(H/3) and W a bs, and the equivalent widths of the H I Balmer 
lines. Colons indicate errors of the order or larger than 40%. 
As explained before, the contribution of the Hg I A4358 line 
flux to the [O ill] A4363 line flux was subtracted assuming 
the average value of the Hg I A4358 line and the number of 
fibers used to get the integrated spectrum of the region. 

The Hq map obtained using 2D spectroscopy allows us 
to determine the total luminosity and the mass of the ionized 
gas within these regions. The extinction-corrected Ha flux 
was determined considering 



/(Ha) = F(Ra) x 10 



c(H/3)x[/(H Q ) + l] _ p/^ a \ x 10 A Ha /2.5 

(2) 

The value of the extinction in the Ha line, Au a , was 
computed in each reg ion considering the value of the 
/(Ha) provided by the ICardelli. Clayton fc Mathisl (|l989l ) 
extinction law and the c(H/3) derived in our analysis, 
Ana = 1.755xc(H/3), and it is tabulated in Table Q] 
The total Ha luminosity was computed assuming a dis- 



tance to IC 10 of 0.8 Mpc, the uncertainty of L(Ha) 
also includes the uncertainty in the distance. We esti- 
mate a total Ha luminosity of (1.18 5±0.064)x 10 36 ergs -1 . 
Using the ICalzetti et al.l (|2007l ) calibration between 
L(Ha) and the star-formation rate (SFR) we derive 
SFR Ka = (6.28 ± 0.34) x 10~ 3 M yr _1 . This value corres- 
ponds to a SFR density of 1.64 x 0.09 M© yr" 1 kpc" 2 , and 
confirms the starburst activity that the center of IC 10 is 
experiencing. 

Finally, the mass of the ionized gas, Mhh, was com- 
uted using the expression provided by IPerez-Monterd 
20021 ). Mhii = 1.485xl0" 35 L(Ha)x(n e /100), where n e is 
the electron density (in cm -3 ) derived from our spectro- 
scopic data (which was always 100 cm -3 ) and the result is 
obtained in solar masses. The derived values of L(Ha) and 
Mhii of each region are also included in Table [T] We es- 
timate a total ionized gas mass of (1.76 ± 0.10) x 10 Mq 
within all the observed field. Table [1] also lists the ionized 
gas mass per surface area, Mhii /area. It has a minimum 
value of 1.10 Mq pc -2 in the hole but it reaches the max- 
imum values in the strong star-forming region, which has 
MHii/area=15.3 M pc -2 . 



4.1 Reddening map 

Left panel of Fig. [4] shows the map of the reddening coef- 
ficient derived in IC 10 using the Ha/H/3 ratio. c(H/3) has 
values between 1.0 and 1.8 dex. These values correspond to 
an extinction in the l/-band, Ay, between 2.15 and 3.86 mag 
or, equivalently, a color excess, E(B — V), between 0.69 
and 1.25 mag. The integrated value of the extinction in the 
V^-band is Av = 3.05 ±0.11 mag. These values nicely match 
previous values of the extinction given in the literature (e.g., 
Lee, Freedman & Madore 1993; Vaccaet al. 2007; Lozinskaya 
et al. 2009). These values confirm the high extinction exist- 
ing in the direction to IC 10, but because of its very low 
Galactic latitude we know that a very important fraction of 
the observed extinction is due to the Galactic ISM. In any 
case, assuming that the Galactic extinction is mainly con- 
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0.8 1.0 1.2 1.4 1.6 
c(H/S) from Ha/Hfi [ dex ] 
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c(H/3) from Ha/H/S assuming Wobs = 2.2 A [ dex 



Figure 5. Comparison between the c(H/3) derived from the Hct/H/3 ratio (x-axis) and that derived from the H7/H/3 ratio (y-axis). The 
x = y is plotted with a discontinuous red line. (Left diagram) c(H/3) derived no considering H I Balmer line absorption. (Right diagram) 
c(H/9) derived assuming an absorption of W a t, s =2.2 A in the H I Balmer lines. 



stant in all the observed region, we appreciate that the dust 
within the star-forming region has a rather inhomogeneous 
distribution. In particular, we note the the lower values of 
c(H/3) are located near the WR stars, at the SE of the bright- 
est H 11 region, that show c(H/3) between 1.0 and 1.3 dex. 
The maximum values of the reddening coefficient are found 
at the NW of the brightest star-forming region, which has 
values between 1.5 and 1.8 dex. Other areas (at the south 
of the hole and in the knot a at the SW corner) also show 
a high extinction. The average value we determine for all 
the observed region is c(H/3) = 1.42±0.05 dex. This peculiar 
distribution suggests that the dust around the massive stars 
has been evacuated by interstellar winds and/or radiation 
pressure or partially destroyed. Although this effect has not 
been investigated in detail, there are indications of decreas- 
ing values of the extinction around young stars in the Orion 
Nebula Cluster (|Hillenbrandll 19971) . 



We compared the reddening coefficient derived from the 
Hq/H/3 ratio with that estimated using the H7/H/3 ratio in 
those fibers in which H7 is properly measured (52% of the 
spaxels). The left panel of Fig. [5] plots the c(H/3) values de- 
rived in each case. As we see, although the correspondence 
between both values is relatively good, the c(H/3) derived 
from the H7/H/? ratio is systematically lower than that com- 
puted using the Hq/H/3 ratio. This effect is consequence of 
the H 1 Balmer line absorption of the underlying stellar pop- 
ulation. The right panel of Fig. [5] plots the c(H/3) values 
derived in each case when considering a H I Balmer line ab- 
sorption of Wai, s =2.2 A, the value we derived iteratively for 
the integrated spectrum of the observed field (see Table [TJ. 
Indeed, the correspondence between both c(H/3)s is much 
better when W a bs is taken into account, although the dis- 
persion of the data is slightly higher because of the higher 
uncertainty in both the H7 fluxes and equivalent widths. 



4.2 Physical conditions of the ionized gas 

Because of the faintness of the auroral lines and the conta- 
mination of the Hg 1 A4358 telluric line, it was not possi- 
ble to derive a temperature map from our PMAS data. We 
used the [S 11] A6716/A6731 ratio to derive a bidimensional 
map of n e , assuming a T e =10500K (the electron temper- 
ature derived for the ionized gas, see below). Lower val- 
ues of the [S 11] A6716/A6731 ratio indicate higher values 
in n e . We used the temden task of the nebular package of 
iraf l|Shaw fc Dufourlll995h with updated atomic data (see 
iGarcfa-Roias et al.ll2005T ). This map is shown in the right 
panel of Fig. U Purple and blue points indicate the places 
where the electron density is in the low limit, n e <100cm -3 
([S 11] A6716/A6731 ratios between 1.40 and 1.32). As we 
see, these colors are the predominant in the observed region. 
However, we note some places with electron densities be- 
tween 100 and 400cm~ 3 ([S 11] A6716/A6731 ratios between 
1.32 and 1.10). Hence, these regions have a higher gas den- 
sity than that found in the surrounding ISM. It is important 
to note that the three regions with the highest electron den- 
sity are precisely located in the places where we find the 
maximum of the reddening coefficient (left panel in Fig. Q. 

For the integrated spectra, we derived the flux of the au- 
roral [O ill] A4363 line by subtracting the expected flux of 
the telluric Hg I A4358 line considering the number of fibers 
used in each region. We then applied the temden task of 
the nebular package of iraf to derive the T e [O 111] in each 
region. For AP 2, which integrates the flux of the faintest 
regions of the nebula, the value derived for the [O ill] A4363 
line is just a residual value. In all cases, T e [0 ill] is around 
10500 K, in agreement with previous resu l ts found in IC 10 
(i.e., lGarnettlll990l ; ICrowther et al.ll2003l ; lLozinskava et all 
2009), so we are quite confident in the validity of our 
decontamination of Hg I emission. Using the relat i on be - 
tween T e [0 III] and T e [0 11] provided by iGarnettl (|l992h . 
T e [0 n]=0.7xT e [O m]+3000, we have estimated the elec- 
tron temperature for the low-ionization potential ions. We 
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[0 III] X5007 /H0 mop (red, corrected) + Ha contours 




He I A5876 / Ha mop (red. corrected) + Ha contours 




SBA [orcsec] 
[N II] A6583 / Ha map (red. corrected) + Ha contours 
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Figure 6. Emission line ratio maps with Ha contours overlaid into the IC 10 PMAS field of: (tot left) [O III] A5007/H/3, (top right) 
He I A5876/Ho, (bottom left) [N n] A6583/Hc?, (bottom right) [S II] AA6716,6731/Ha. All maps have been corrected for reddening using 
the c(H/3) value derived from the Hcr/H/3 ratio. Symbols are the same that in Fig. [2] 



compile the electron temperatures derived in each of the 
regions in Table [2] This table also lists the electron density 
computed for each region using the [S n] A6716/A6731 ratio. 
In all cases, this ratio provided a density which was below 
the low-density limit, n e <100 cm -3 , and hence we always 
considered n e = 100 cm -3 . 



4.3 Nature of the ionization 



ob- 



Figure [6] shows the emission line ratio maps we 
tained for the analyzed region in IC 10. All maps 
have been corrected for reddening using the c(H/3) value 
derived from the Ha/H/3 ratio. We show the maps 
corresponding to [O m] A5007/H/3, [N n] A6583/Ha, 
([S n] A6716+A6731)/Hq, and He I A5876/Ha. White colors 
indicate regions in which the studied emission line was not 
detected or it was too noisy. 

As we see, the maps are showing the ionization struc- 
ture of the giant H II regions in the center of IC 10. 



That is specially evident in the brightest star-forming re- 
gion [HL90] 111c and in the H II region located at the SE, 
[HL90] Hid, both show the maxima of the [O ill] A5007/H y S 
ratio at their centers, while the maxima of the low- 
ionization potential ion ratios is located surrounding those 
central bright areas. This effect is clearly noticed in the 
[O ill] A5007 / [O II] A3727 ratio map plotted in FigH which 
shows the higher values at the center of the strong star- 
forming regions. This result is very important and should be 
taken into account, as we will see below, when using empiri- 
cal calibrations to get the chemical abundance of the ionized 
gas. It is also interesting to note that the highest values of 
the [O ill] A5007/H/3 ratio is close to the maximum of the 
Ha emission, but they are not completely coincident. 

We have also analyzed the diagnostic diagrams plot- 
ting two different excitation line ratios for classifying the 
excitation mechanism of the ioni z ed ga s, as initially pro- 
posed by IVeilleux fc Osterbrockl (|1987m . Figure [8] plots 
the typical [O ill] A5007/H,S versus [N n] A6583/Ha and 
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X5007 / [0 111 A3727 (red. corrected) + Ha contours 
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Figure 7. Map of the dereddened [O III] A5007 / [O II] A3727 ra- 
tio computed in the intense star-forming region [HL90] 111 within 
the starburst galaxy IC 10. The contours of the Ho emission are 
overlaid. Symbols are the same that in Fig. [2] 



[O in] A5007/H/3 versus ([S n] A6716+A6731)/Ha diagrams 
for all the spaxels for which we have a good measurement 
of these lines. In these diagnostic diagrams, H n regions and 
starburst galaxies lie within a narrow band, but when the 
gas is excited by shocks, accretion disks, or cooling flows 
(in the case of AGNs or LINERs) its position is away from 
the locus of H n regions. We used the anal ytic relations 
given by the photoionized models provided by iDopita et al.l 
(2000) for extragalactic H n regions (that assume instanta- 
neous star-formation within star-forming regions) and the 
iKewlev et al.l {2001) models for starburst galaxies (which 
consider continuous star formation and more realistic as- 
sumptions about the physics of starburst galaxies). In both 
cases, photoionization models with very different metallic- 
ities (between 0.1 and 3 Zq) and ionization parameters 
lie always under these dividin g lines . Indeed, the dividing 
line given by the IKewlev et aD (|200lh models represents an 
upper envelope of positions of star-forming galaxies. The 
left panel of Fig. [5] includes the empirical relation between 
th e [O hi] A5007/H / 3 and the [N n] A6583/Hq provided 
by iKauffmann et all l|2003h analysing a large data sample 
of star- forming galaxies from the Sloan Digital Sky Survey 
(SP SS; [York et al.l l 20QPj) . 

IKewlev et al.l (|2006h suggested th at those obje c ts be - 



2001 



2003) 



tween the theoretical line compute d by Kewlev et alj 
and the empirical line found by IKauffmann et ak 
may be ionized by both massiv e star s and shocks. How- 
ever, iPerez-Montero fc Continil (|2009h showed that ob- 
jects located in this area may also be pure star-forming 
galaxies but with a high N content. In our case, all 
data points (even those spaxels for which a possible 
N enrichment has been detected, see § 5.1) lie within 
this area, but below the IKewlev et al.l {2001) theoreti- 
cal line. This clearly indicates that photoionization is the 
main excitation mechanism of the gas and that there 
is very little evidence for a significant contribution from 
shock excitation. The position of the data points in the 



o 0.0 
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Figure 8. Comparison of some observational flux ratios obtained 
for all available spaxels with the diagnostic diagrams proposed by 
Dopita et al. (2000), blue continuous line (D00), and Kewley et 
al. (2001), red discontinuous line (K01). The left panel also shows 
the empirical relation provided by Kauffmann et al. (2003) with 
a dotted-dashed dark yellow line (Ka03). 



[O in] A5007/H,S versus [N n] A6583/Ha diag ram perfectly 
match es with the predictions given by the IDopita et alj 
(2000) models. In the case of the [O ill] A5007/H/3 ver- 
sus ([S n] A6716+A6731)/Hq diagram, the data points also 
agree well with these models, but the data te nd to move to 
the pr edictions of the upper envelope given by IKewlev et al.l 
<|200lf > in points with log([0 in] A5007/H^)<0.5 and 
log([S n] A6716+A6731)/Ha> -0.8. 



4.4 Chemical Abundances 

Because of the lacking of an individual estimation of the 
electron temperature, we are not able to map the chemical 
abundances using the direct method. We will discuss in i]5.2l 
the effects of using empirical calibrations to map the che- 
mical abundances of an H II region or a starburst galaxy 
using 2D spectroscopy. Here we analyze the chemical abun- 
dances obtained from the integrated spectra of the areas we 
described in Fig. [2] 

We used the iraf package nebular to derive ionic 
abundances of N+, + , 2+ , S + , Ne+ 2 , and Ar +2 from the 
intensity of collisionally excited lines. We assumed a two- 
zone scheme for deriving the ionic abundances, adopting 
T e [0 ill] for the high ionization potential ions 2+ , Ne+ 2 , 
and Ar +2 , and T e [0 n] for the low ionization potential ions 
+ and S + . The errors in the ionic abundances have been 
calculated as a quadratic sum of the independent contribu- 
tions of errors in flux, n e and T e . 

To derive He + /H + , we used three observed lines of He I 
at AAA4471, 5876, and 6678, weighted by 1:3:1. Case B emis- 
sivities for the T e [0 ill] and n e assumed for each region, 
were taken from the colli s ion-le ss (low-density limit) calcu- 
lations by iBauman et 3 (|2005l) using an on-line available 



10 A.R. Lopez-Sanchez et al. 



Table 2. Physical conditions and chemical abundances of the ionized gas for the regions analyzed in IC 10 using our PMAS data. Three 
dots indicate that the abundance could not be derived because of the non-detection of some emission lines. 



Region 


[HL90] 111 


API (c) 


AP2 (hole) 


AP3 (d) 


AP4 (WR 24) 


Fiber 90 


Fiber 92 


rp t|t] \T/*i 

i e [U III] [rvj 


lUOUUztOUU 


lUOUUztoUU 




lU/UUztoUU 


1 nc^nn-Li nnn 
lUOUUztlUUU 


1 nf;/in4- 1 nnn 
lUOUUztlUUU 


1 nt^nn-Ucnn 


T„ ro III fKl 


10400±400 


10400±400 


10400±700 


10500±600 


10400±700 


10400±700 


10400±600 


Yi „ [nn ^1 


100 


100 


100 


100 


100 


100 


100 


1 0_Llnrr/'n+ /'P'+ A 


7 Q^-Ln 1 


/ . / o±U. 1U 


7 fio^n oa 
t .oy±U.Z4 


7 C/1 _i_n 1 1 


7 7/1 4~n 1 o 


7 7n^n on 


7 RK-i-Cl 1 Q 
/ .Du±U. lo 


12+loe:('0 ++ /H + ) 


7.96±0.06 


8.03±0.07 


7.94±0.10 


8.01±0.07 


8.02±0.08 


8.05±0.09 


8.07±0.08 


12+log(0/H) 


8.26±0.09 


8.22±0.08 


8.22±0.17 


8.23±0.09 


8.20±0.12 


8.21±0.13 


8.21±0.09 


log (O++/O+) 


0.00±0.13 


0.25±0.12 


0.04±0.25 


0.17±0.14 


0.28±0.20 


0.36±0.21 


0.42±0.15 


12+log(N+/H+) 


6.54±0.05 


6.33±0.05 


6.49±0.15 


6.43±0.05 


6.49±0.10 


6.24±0.11 


6.19±0.06 


12+log(N/H) 


6.84±0.08 


6.78±0.09 


6.81±0.20 


6.90±0.10 


6.95±0.17 


6.75±0.18 


6.76±0.12 


log(N/0) 


-1.42±0.12 


-1.45±0.11 


-1.41±0.24 


-1.41±0.13 


-1.26±0.20 


-1.46±0.21 


-1.45±0.15 


12+log(S+/H+) 


5.77±0.04 


5.50±0.04 


5.75±0.09 


5.67±0.05 


5.54±0.09 


5.41±0.09 


5.40±0.05 


log(S+/0+) 


-2.19±0.16 


-2.28±0.14 


-2.14±0.33 


-2.17±0.17 


-2.20±0.28 


-2.29±0.29 


-2.25±0.18 


log(N+/S+) 


0.77±0.09 


0.83±0.09 


0.74±0.24 


0.76±0.10 


0.95±0.19 


0.83±0.20 


0.79±0.11 


12+log(Ne++/H+) 


7.50±0.16 


7.45±0.19 




7.50±0.25 






7.38±0.21 


12+log(Ne/H) 


7.80±0.16 


7.65±0.18 




7.73±0.25 






7.52±0.21 


log(Ne/0) 


-0.46±0.23 


-0.58±0.24 




-0.51±0.30 






-0.69±0.27 


12+log(Ar++/H+) 


5.90±0.08 


5.92±0.09 


5.81±0.16 


5.94±0.10 


5.90±0.13 


5.90±0.13 


5.96±0.09 


log(Ar++/0++) 


-2.06±0.14 


-2.11±0.16 


-2.13±0.26 


-2.07±0.17 


-2.12±0.21 


-2.15±0.22 


-2.11±0.17 


12+log(He+/H+) 


10.99±0.04 


10.98±0.03 


11.03±0.11 


10.99±0.04 


11.09±0.05 


10.97±0.09 


10.99±0.03 



codfl The collisional to recombination contribution was es- 
timated from iKingdon fc Ferlandl lll995lK using the inter- 
polation formula provided by iPorter et al.l l|2007l ). The ef- 
fective recombinat i on co efficients for H + were taken from 
IStorev fc Hummer! (|l995l ). 

We computed the total abundances of O and N, adopt- 
ing O/H = 0+/H+ + +2 /H+ to determine the total oxy- 
gen abun dance and the standard ioniz ation correction factor 
(ICF) bv lPeimbert fc Costerol (|l969h . N/O = N+ /0+, to 
compute the total nitrogen abundance. For those regions in 
which we have a detection of the [Ne ill] lines, we computed 
the total neon abundance applying the ICF proposed by 
iPeimbert fc Costerol (|l969T ). Ne/O = Ne+ 2 /0+ 2 . 

We derive a total O abundance of 
12+log(O/H)=8.26±0.09 for all the observed region. 
The derived N/O and Ne/O ratios, log(N/0)=-1.42 ± 0.12 
and log(Ne/0)~ —0.42, are the typic al values expected for a 
galaxy with the metallicity of IC 10 (llzotov fc Thuanl l 19991 : 
llzotov et af]|2004l ; lLopez-Sanchez fc EstebanlbOlObl ). These 
results are very similar in the regions AP 1 ([HL90] 111c, 
the brightest star- forming region), AP 2 (the Ha hole) and 
AP 3 QHL90] 11 Id), but they slightly differ in the AP 4 
region, which only considers the emission in those two 
fibers where the WR stars are detected. We will discuss this 
intriguing result in Sect 15. ll 

4.5 Wolf-Rayet content 

We have clearly detected the blue WR bump in two adja- 
cent fibers (spaxels [11,8] and [11,9]). This position coincides 
with one of the WR stars catalogued within [HL90] 111c 
(see FigQ] right) in the most r ecent WR stars catalogue 
of IC 10 (|Crowther et alj I2003T ), [MAC92] 24C. We have 

2 Available at http://www.pauKy.edu/~rporter/j-resolved 



also a tentative detection of the nebular He II A4686 line in 
spaxel [11,10], which corresponds to the position of the WR 
[MAC92] 24B. However, we do not observe any WR feature 
at the position of [MAC92] 24A (spaxel [10,10]) or in the 
position of RSMV 23 (spaxels [2,14] and [2,15]), although 
the S /N ratio of this last region is relatively low. 

We have added the spectra of the two spaxels where the 
blue WR bump is observed (AP 4) and carefully analysed 
the properties of the integrated spectrum. The emission line 
fluxes are listed in Table [1] and the derived physical and 
chemical properties of the ionized gas are given in Table [2] 
For comparison, we have also analyzed two nearby fibers, 
the adjacent position to the detection of the blue WR bump 
(F90, which corresponds to the spaxel [11,10], where there 
is a tentative detection of the nebular He II A4686 line), and 
the fiber with the maximum of the Ha emission (F92, which 
corresponds to the spaxel [11,12]). 

In our case, the blue WR bump we detect in AP 4 is 
essentially constituted by the stellar, broad, He II A4686 
emission feature. Figure [9] shows a detail of this spec- 
trum in the 4620-4750 A range, showing a prominent 
bump. We do not detect the C m/C IV A4650 or the 
C iv A5808 (the red WR bump) an d hence we agree with 
the results found by other aut hors (|Crowther et alj 120031 ; 
IVacca. Sheehv fc Grahamll2007f l that the observed WR fea- 
tures are originated by a WN star. Both the strength of the 
broad He II A4686 emission and the absence of the N v A4604 
emission indicate that this is a late- type WN star (WNL). 

We have followed the method described in 
lLopez-Sanchez fc Estebanl l|2010al ) to estimate the flux 
of the blue WR Bump. We fitted a broad and a nar- 
row Gaussian for the stellar and nebular He II A4686 
lines. Figure [9] shows the result of this fitting, in- 
cluding the residual spectrum after subtracting the 
best-fit model to the observed spectrum. The best- 
fit model values are _Fl4 y -H"M H oiibroad=8.12±0.76 A, 
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Figure 9. (Top) Spectrum, in the 4620—4750 A range, obtained 
adding the two fibers (spaxels [11,8] and [11,9]) in which the blue 
WR bump is observed (AP 4). The red dotted line represents 
the position of the He II A4686 emission line, the black dashed 
line represents the continuum fit. The best fit to the observed 
spectrum (blue continuous line) is shown with a black continuous 
line, and it is composed by the superposition of a broad and 
a narrow He II A4686 lines, following the method explained in 
Lopez-Sanchez & Esteban (2010a). (Bottom) Residual spectrum 
after subtracting the best-fit model to the observed spectrum. 



fHeiibroad = (8.76±0.79)xl(r 16 ergs -1 cm -2 for the 
broad component and FW-ffMHciinarrow=5.4±0.9 A, 
JHeiinarrow = (6.1±1.0) x 10 -17 ergs -1 cm -2 for the narrow 
component. 

Using the reddening coefficient derived for AP 4 
and its relation with the extinct i on at A4686 [assum- 
ing the ICardelli. Clayton fc Mathisl l| 19891 ) extinction law 
with 7?v=3.1], Aa4686 = 2.625xc(H/3), and considering 
the distance to IC 10, we derive a total luminosity 
of L HeIIMoad = (1.35 ± 0.17) x 10 36 ergs -1 for the broad 
He II A4686 emission. We can use this value to compute 
the total number of WR stars in t h e burs t. For this, and 
following lLopez-Sanchez fe Esteban (2010a), we consider a 
metallicity-dependence of the WR luminosities. Using the 
relation between the luminosity of the He il A4686 line 
and the oxygen abundance (their Eq. 7), we estimate that 
a single WNL with the metallicity found in IC 10 [which 
we assume to be 12+log(0/H)=8.25], has a luminosity 
of L W NL(He n A46 86)= 1.27 x 10 36 ergs" 1 . Hence, the 
blue WR bump observed in this position of IC 10 is con- 
sistent with be ing produced by a single WNL star and 
not by four, as IVacca. Sheehv fc Grahaml l|2007l )suggested 
from their optical/NIR images. This conclusion is es- 
sentially the same in the case of considering the lumi- 
nosity of a WNL star with a solar metallicity to the 
broad A4686 flux, £yy NL (H e n A4686)= 1.6 x 10 36 ergs" 1 
(ICrowther fc Hadfieldl i2006). which gives 1.19 WNL stars 
within this starburst. 

To estimate the WR/(W R+Q) ratio we also fol- 
low th e prescriptions given by lLopez-Sanchez fc Estebar] 
l|2010al ). We determine this ratio within the intense 



star-forming region. The total number of ionizing pho- 
tons is Qtot = Lh/3/4.76 x 10" 13 = 2.84 x 10 50 s" 1 , an d the 
r?(t)=Q7V/Q parameter (|Vacca fc Contilll992l ; IVaccalll994l ; 
ISchaerer fc Vaccal 1 19981 ) which depends on the age of 
the burst, is t] ~ 1 . Hence , applying Eq 10 of 
lLopez-Sanchez fc Estebar] |2010al ) considering Nwnl~1 
and Nwce=0, we estimate a total number of O stars of 
26 and a WR/(WR+0) ratio of 0.037. This value agrees 
with that determined in other WR galaxies with a similar 
oxygen abundance (see Fig. 5 in lLopez-Sanchez fc Esteban! 
l2010al ). 

Fiber 90 (spaxel [11,10]) does not only shows a tentative 
detection of the nebular He II A4686 line, but it also shows 
a clear broad Ha component underlying the narrow, neb- 
ular Ha emission. These kinds of featur es are rather com- 
mon in extragalactic H II region s (e.g.. [Firpo et al.l [ 2005 . 
2010) and starburst galaxies (e.g.. iMendez fc Estebanll 19971 : 
Homeier fc Gall agher 1999|). Recent detailed analyses of 
nearby starbursts us i ng IFU spectrosc opy [e.g., NGC 1569, 
IWestmoauette et al.l (l2007al); M 82 , IWestmocmette et all 



I 2007bl ); Mkn 996. 1 James et all (|2009r i] are also revealing the 
existence of broad components underlying the bright nebular 
lines that are originated in turbulent movements and shocks 
of the ionized gas. Our analysis of the emission-line profile s 
within the center of NGC 5253 (|L6pez-Sanchez et al]|2007f ) 
also indicated that they were composed by a narrow an d 
a broad components (see also iMonreal-Ibero et al1l2010l ). 
Furthermore, the analysis of their chemical abundances re- 
vealed that the broad component of the knot B of NGC 5253 
seems to contain the localized N enrichment in this area of 
the galaxy (see their §9 and their Table 10). Very different 
abundance results between indi vidual line compon ents were 
recently reported in Mkn 996 (jjames et al.l 2009). Because 
of the lacking of spectral resolution, we cannot do this analy- 
sis here. However, we have followed the same procedure used 
before to fit the nebular and stellar He II A4686 lines to fit the 
narrow and broad components of the Ha line in this fiber. 
We estimate that the extinction-corrected flux contained in 
the broad Ha component is (8.3 ±0.7) x 10~ 15 ergs" 1 cm -2 , 
which is ~11 % of the flux measured in the narrow line. This 
value indicates that this component contains a mass of about 
~9.5 M e . The estimated FWHM of the broad Ha compo- 
nent is ~40 A, which indicates a velocity of ~1830 kms -1 . 
This high value suggests that the flow should be related to a 
remnant of a recent SN explosion that has still not produced 
an observable chemical enrichment. 



4.6 Age of the starburst 

Figure \W\ shows the map of the equivalent width of the Ha 
emission, Who, derived in the intense star-forming region 
[HL90] 111 from our PMAS data. As we see, the maximum 
values of Who, are found at the west of the central star- 
forming region, which shows Who, between 850 and 1000 A. 
Using the STARBURST 99 stellar population synthesis 
model (|Leitherer et al.l |l~999) with an instantaneous burst 
with a Salpeter IMF slope (a=2.35) with M up =100 M© and 
Z=0.008 Zq [which corresponds to an oxygen abundance of 
12±log(0/H)=8.26 considering a solar ox ygen abundance of 
12±log(0/H)=8.66 (|Asplund et al.ll2005l ) we derive an age 
between 3.3 and 4.2 Myr for the most recent starburst. The 
eastern areas of this intense star-forming region shows val- 
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Figure 10. Map of the equivalent width of the Ha emis- 
sion, Whqi derived in the intense star-forming region [HL90] 111 
jHodge fc LeeHl990t) from our PMAS data. The contours of the 
Ha emission are overlaid. Symbols are the same that in Fig. [2] 
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Figure 11. [N n] A6583 / ( [S n] A6717 + [S n] A6731) line 
intensity ratio map derived for the strong star-forming region 
[HL90] 111 within the starburst galaxy IC 10. The contours of the 
Ha emission are overlaid. Symbols are the same that in Fig. [2] 



ues around Wh«~600-700 A, which are translated to an 
age of ~4.7 Myr. Similar values are found at the central re- 
gion of the smaller star-forming region at the east (knot d). 
The arc-like structure at the south has Whq~300-500 A. Al- 
though the presence of Ha absorption make underestimate 
the actual equivalent widths of the emission lines (and hence 
the age of the last star-formation event is overestimated) we 
consider that in our case their effect in Wu a is very small 
because of the strength of this emission line. Indeed, the 
average value of W a bs derived here (see §4.1) is only 2.2 A. 
However, as we discussed before, the effect of the H I line ab- 
sorption should be considered when studying the rest of the 
Balmer lines and, particularly, to derive a more appropriate 
value of the reddening coefficient. 

It is interesting to note that the places in which the 
WR features are detected have relatively low values of Wn a , 
^100 A, which yields to an age of ~6.2 Myr for the most 
recent star-formation event. That age is older that that we 
should expect for a WR star, as the most massive, luminous 
and hot O stars evolve to the WR phase between 2 and 5 Myr 
since their birth, spending only some fe w 10 5 years before 
they explode as type Ib/Ic supernovae (|Mevnet fc Maederl 
2005). Hence, the star-formation activity may not be strictly 
instantaneous but it has probably had some continuous rate. 
Maybe, the detected WR star was formed somewhat later 
after the majority of the stellar population in the area. Fur- 
thermore, as we are considering very small spatial scales 
(1"— 3" ~3.9-11.6 pc), the WR star may currently be some 
distance from its birth plac^f). 



5 DISCUSSION 

5.1 A possible localized chemical pollution 

3 A star with a proper motion of 20 km s~ 1 will move ~10 pc in 
1 Myr. 



by a Wolf-Rayet star? 

The only confirmed case of localized chemical pollution by 
WR stars in a starb urst region to date i s the center of the 
BCD G NGC 5253 llWaTsh fc Rovl [l989l; iKobulnickv et all 
Il997t lLopez-Sanchez et al.ll2007r i." This object shows a clear 
enhancement in the N/H ratio (and also in the He/H ratio) 
in a small region, which is consistent with the chemical pollu- 
tion expected by the contribution of just a few (1 -2) evolved 
(WR) massive stars (|L6pez-Sanchez et alj|2007h . There are 
some studies that also suggest that starbursts in which WR 



stars are detected may show an enhancement in their N/O 
ratio llPustilnik et al. [120041; I Brinchmann. Kunth fc Durret l 
120081; I James et al.ll2009l ; lLopez-Sanchez fc Estebanll2010bl) 



i^en t 2D spectro sTopTc al data |Pere^Montero^M ) 
also report very high N/O ratios in BCDGs with WR stars. 
However, in these cases the chemical pollution by the WR 
stars detected in the central region seems not to be able to 
explain the large enhanc ement of nitrogen found throughout 
the B CDG. Similarly, |X morin. Perez-Montero fc Vflchezl 
(2010) used SDSS data to find that, at a given metallicity 
the so-called green pea galaxies display systematically large 
N/O ratios compared to normal galaxies, but these spectra 
are not deep enough to detect WR features. 

In this case, our 2D spectroscopic PMAS maps of the 
strong star- forming region [HL90] 111 within IC 10 sug- 
gests the detection of a localized high N/O ratio at the po- 
sitions where the WR features are detected (AP 4). Fig- 
ure [IT] shows the [N n] A6583 / ( [S n] A6717 + [S n] 
A6731) line intensity ratio map. Because of the faintness 
of the [O II] A3727 doublet, the [N n] A6583 / [O n] A3727 
line intensity ratio map shows a very low S/N in a large 
number of spaxels where the surface brightness of the neb- 
ular gas is low, and hence it is useless for our purposes. 
Considering only the spaxels with higher S/N ratio, the 
average value of the [N n] A6583 / [O n] A3727 ratio is 
~0.07. Spaxel [11,8], labelled with a cross in all maps and 
where some He II A4686 emission is observed, has a very 
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uncertain detection of the [O n] A3727 doublet. However, 
spaxel [11,9] (labelled with a star and where the highest 
He II A4686 emission is observed) seems to have a slightly 
higher [N n] A6583 / [O II] A3727 ratio, ~0.10, but because 
of the high uncertainty this value is not significant. On the 
other hand, the [N n] A6583 / ( [S n] A6717 + [S n] A6731) 
ratio, which is well measured in almost all spaxels, 
indicates a rather clear increase of its average value 
(~0.85) precisely at the two spaxels where WR features 
are detected. Indeed, spaxels [11,8] and [11,9] show a 
[N n] A6583 / ([S n] A6717+[S II] A6731) ratio of 1.11 and 
1.01, respectively, with a typical uncertainty of ±0.12. 

It is interesting to note that, the He I A5875/Ha ra- 
tio map (top right panel in Fig. [5]) also shows an enhance- 
ment of the emission at the same position. The average He I 
A5875/Ha ratio is ~0.046 but spaxels [11,8] and [11,9] have 
a ratio of 0.071 and 0.058, respectively, with an error of only 
±0.007. It is possible that we are detecting at this position a 
localized enhancement of both nitrogen and helium, which, 
if true, is very probably connected with the detection of WR 
features in those same spaxels. 

We further investigate this very intriguing result via 
the careful analysis of the chemical properties of the inte- 
grated spectrum obtained adding both spaxels (AP 4). As 
control, we also study the adjacent spaxel [11,10] (Fiber 90), 
which has enough S/N ratio to perform such analysis. Ta- 
ble [2] compiles all the chemical abundances and ratios de- 
rived for these two regions. As we see AP 4 possesses, sys- 
tematically, higher N/O, N+/S+ and He+/H+ ratios (-1.26, 
0.95 and 11.09, respectively) than those values derived for 
any other aperture or spaxel, which have average values of 
— 1.44, 0.79, and 10.99, respectively. We also notice that the 
S + /0 + (two elements of similar nucleosynthetic origin) ratio 
remains practically the same for all regions including AP 4, 
showing a average value of —2.22. Hence, although we can- 
not definitively confirm that there is a localized chemical 
pollution in nitrogen in those two spaxels because of the 
relatively high uncertainties, all evidence points towards it. 
We must emphasize that the helium enrichment is indeed 
confirmed within the errors, because the average value of 
12±log(He+/H+) is 10.99±0.03 while that measured in AP 4 
is 11.09 ±0.05. 

With the data collected in tables [T] and [2] we can make 
a rough estimation of the mass of newly created He and N 
(stellar yield) necessary to produce the observed overabun- 
dances in region AP 4. These yields have been computed 
considering the He and N abundances of region AP 1 as 
reference. Therefore, we have assumed that the initial abun- 
dances of the ionized gas at AP 4 are those measured at 
AP 1 and that AP 4 has suffered a localized and very recent 
increment of the He/H and N/H ratios. Finally, the em- 
pirical stellar yields have been obtained multiplying those 
He/H and N/H increments by the mass of ionized gas de- 
rived for AP 4. The computed stellar yields are 28±16 Mq 
and 0.013±0.013 M for He and N, respectively. These 
numbers can be compared with those given in table 12 of 
lLopez-Sanchez et al.l l|2007h . that includes the yields deter- 
mined for regions A and B of NGC 5253 as well as theoretical 
ones from stellar evolution models and empirical determina- 
tions for ring nebulae around Galactic WR stars. The yields 
found in [HL90] 111c (although uncertain) are roughly con- 
sistent with the scenario of chemical pollution by the ejecta 



of a very small number of WR stars. The amount of He ap- 
pears to be rather high, in fact the ratio of He and N stellar 
yields is about 2150, ten times higher t han in the objects in- 
cluded in the aforementioned table of lLopez-Sanchez et al.l 
|2007ft . 

We must say that it is not the first time that a re- 
gion with a probabl y high N/O is reported in IC 10. 
iMagrini fc Goncalves] |2009l ') used Gemini North Multi- 
Object Spectrograph to map the H II regions and Plane- 
tary Nebulae (PN) in the central 5.5' x 5.5' of IC 10. Their 
data suggested that the H II region [HL90] 120, with an 
oxygen abundance of 12+log(0/H)~8.18, has a N/O ratio 
of ~ —0.48 dex, but the other H II regions (which have an 
oxygen abundance between 8.0 and 8.7) show N/O ratios of 
~ —1.25. However, their data have high uncertainties (the 
majority of N/O values agree within the errors, see their 
Fig. 2), and hence this result should be also confirmed. Fur- 
thermore, all the N/O ratios these authors derived for their 
analyzed H II regions within IC 10 seem to be higher than 
those expected for their oxygen abundance. 

As IMagrini fc Goncalves! (|2009h pointed out, the H n 
region [HL90] 120 contains four known WR stars within 20" 
from its center, and hence the origin of the tentative high 
N/O ratio may be some of these WR stars. In our case, 
we notice that our possible detection of a high N/O and 
He/H ratios are found exactly at the same position where 
the WR features are located, indicating that, if real, this 
chemical pollution is indeed very localized (~ 2" ~ 7.8 pc) 
and therefore it should be very difficult to detect in fainter 
and distant galaxies. 

Some studies affirm that the more or less low and 
constant values of the N/O ratio reported in BCGs is 
a consequence of the fact that anything ejected by mas- 
sive stars (whether oxygen by type II supernovae or nitro- 
gen in the winds of WR stars) is too hot to be immedi- 
ately incorporated to the ISM, and it needs to cool before 
mixing with the existi ng nebular gas and app ears in the 
optic a l spectrum (e.g., [K obumicky & Skillman 1996, [l997l. 
1 19981 ; iKobulnickv et all Il997l; iHenrv. Edmunds fc Koppen 



2000; Pilvugin. Thuan fc Vi'chez 



2006J; IMagrini fc Goncalves! 2009 



20031; 



. .van Zee fc Havnes 
Lagos et all 120091'). How^ 



ever, as previously suggeste d bv lPustilnik et al.l (|2004r i and 
lLopez-Sanchez et al.l l|2007h , the presence of localized chem- 
ical enrichment in very young starbursts indicates that 
the timescale of the process should be very short. This 
agrees with the enric hment pattern obser ved in Galac- 
tic WR ring nebulae jEsteban et alj ll992T). The s tatisti- 
cal analysis performed by iBrinchmann et ail ((2008) com- 
paring WR and non-WR galaxy data from SDSS observa- 
tions also yields to a similar conclusion. The recent, more 
det ailed, analysis of a sample of WR galaxies performed 
by lLopez-Sanchez fc Estebanl |2010bl ) also reported high 
N/O ratios in 6 WR galaxies. These authors indicated that 
the high N/O values should dilute in a quick period be- 
cause of both the decreasing of the star-formation activity 
and the releasing of new, fresh material which is dispersed 
and mixed with the existing gas of the galaxy. Further- 
more, the effects of the WR star s into the ISM are also di- 
luted because of aperture effects (|L6pez-Sanchez fc Estebanl 
2009T). Only very careful analys es of such systems (e.g.. 



Lopez-Sanchez fc EstebanlfcoiOrl ) will help to solve the puz- 
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zle of the high N/O ratios sometimes found in starburst 
galaxies. 

5.2 Use of empirical calibrations in 2D 
spectroscopy data 

We also study the results provided by the so-called 
empi-rical calibr at ions of the oxygen abundance [see 
IKewlev fc Ellison! (|2008h and lLopez-Sanchez fe Estebanl 
(2010bJ) for recent reviews] when analyzing data provided 
by 2D optical spectroscopy, which in many occasions have 
a good determination of all important nebular lines in- 
volved in the most common empirical calibrations (i.e., 
[O n] A3727, H/3, [O III] AA4959.5007, Ha, [N n] A6583). 
Ratios between these bright emission lines define the pa- 
rameters -R2 3, P, y, N2, and O3N2 [se e defini tions in Ap- 



pendix A of Lopez-Sanchez fc Estebanl l|2010bT )]. Empirical 
calibrations are very commonly used to get a rough estima- 
tion of the metallicity of star-forming galaxies when the elec- 
tron temperature of the ionized gas cannot be determined 
(i.e., there are no detections of the faint auroral lines). 

Although some precautions must be taken into ac- 
count when analyzing star-forming galaxies w ith inter- 
medi ate metallicities, 7.9<12+log(0/H)<8.4, (|Yin et all 
2007), it seems that the combination of the R23 and 
an excitation parameter (P or y) have the necessary 
information for the determination of accurate abun- 
dances in extragalactic H 11 regions. The comparison 
between the oxygen abundances provided by empiri- 
cal calibrations and those determined via th e dire ct 
method performed by lLopez-Sanchez fc Estebanl (2010 
concluded that the Pi lyugin-method |PilvuginTl2001a 



IPilvurin fc Thuanll2005h is nowadays the best suitable em- 
pirical calibration to derive the oxygen abundance of star- 
forming galaxies. However, the results provided by empiri- 
cal calibrations based on photoionization models ( McGaughl 
Il99ll : iKewlev fc Dopital 120021 ; iKobulnickv fc Kewlevl I200I ) 
are systematically 0.2-0.3 dex higher than those found fol- 
lowing th e direct method. This result has been later con- 
firmed bv lMoustakas et al.l (|2010l ). who found differences of 
~0.6 dex betwee n both kinds of method. On the other hand, 
lYin et all |2007l ) indicated that the N 2 and the O3N2 in- 
dices are useful for calibrating metallicities of galaxies with 
12+logfO/H)<8.5. lLopez-Sanchez fc Estebanl feoiObl ) con- 
cluded that the relations between t hese two parameter and 
the oxygen abundance provided by iPettini fc Pagel (|2004l ) 
give acceptable results for objec ts with 12+log(O/H)>8.0. 
IPerez-Montero fc Continil 120091 ) established that empirical 
calibrations using the N2 parameter are not valid for objects 
with 12+log(O/H)<8.0. However, the effects of such empir- 
ical calibrations using 2D spectroscopic data have not been 
explored so far. 

Figure [12] shows the oxygen abundance maps of the 
observed star- forming region [HL90] 111 within IC 10 
derived f rom four of the most common empirical cali- 
brations: IPilvuginl d2001a t]), which involv e the R23 and 
the P parameters; IKobulnickv fc Kewlevl (|2004T l. which 
is based on the phot oionization models presented by 
IKewlev fc Dopital j2002) and considers the R23 and the y 
parameters to iteratively derive the oxygen abundance and 
the ionization parameter q; the bicubic relation between 
the oxygen abundance and the N2 parameter derived by 



IPettini fc Page! l|2004l ); and the linear fit between the oxy- 
gen abundance and the O 3N2 parameter also determined 
by IPettini fc Pagel (2004). As the oxygen abundance de- 
termined using the direct method in IC 10 is found be- 
tween the low-metallicity [12+log(0/H)5C8.1] and the high- 
metallicity [12+log(Q/H) ^8.4] branches, w e assumed the 
average value between the IPilvuginl (|2001al ibh calibrations, 
but alw ays the metallicity derived for the low branch follow- 
ing the IKobulnickv fc Kewlevl (|2004h calibration. The four 
panels shown in Figure [12] have the same scale in the oxygen 
abundance, between 8.0 and 8.4. 

Th e oxygen abundance map provided by the IPilvuginl 
(2001a. b) calibration (top left panel in Fig. ll2|) shows a more 
or less constant value, although it clearly has lower values 
in the high-ionization regions [12+log(0/H)~8.13] that in 
the low-ionization regions [12+log(0/H)~8.21], but both 
agree within the errors (±0.10 dex) and, although slightly 
lower, with the value determined using the direct method 
[12+lo g(O/H)=8.26±0.09l. If using the IPilvugin fc Thuanl 
l|2005l ) calibration, the resulting metallicity map has a sim- 
ilar structure to that shown in top left panel of Fig. [T^] but 
the derived values are systematically ~ 0. 05-0. 10 dex low er 
than those determined following the IPilvuginl (|2001al lbh 
method. In any case, as we will discuss below it is evident 
that there is some effect of the ionization degree on the deter- 
mination of the oxygen abundance and that this calibration 
seems to be not suitable for looking for chemical differences 
smaller than ~0. 10-0. 15 dex using 2D spectroscopic maps. 

The situation is somew hat similar when using the 
IKobulnickv fc Kewlevl (|2004l ) calibration (top right panel in 
Fig. I12[) . As we see, this calibration provides oxygen abun- 
dances between 8.04 and 8.38, with the lower values lo- 
cated again in the high-ionization regions. Spaxels which 
have higher oxygen abundances following this method are 
typically located in the external areas of the star-forming 
regions, where the ionization degree is lower and the uncer- 
tainty in the [O 11] A3727 flux is high er. In any case, we no- 
tice th at the results derived from the IKobulnickv fc Kewlevl 
(2004) calibration have a higher dispersion that those ob- 
tained following the Pilyugin method, and therefore this 
calibration should be used with caution when analyzing the 
chemical abundances of 2D spectroscopic maps. 

The ionization structure of the giant H 11 regions we 
are observing with our PMAS data is indeed playing an im- 
portant role in the determination of the oxygen abundances 
using empirical calibrations. That is very evident when us- 
ing the N2 and O3N2 parameters, which strongly depend 
on the ionization degre e of the gas. The oxyge n abundance 
map obtained using the IPettini fc Pagell (|2004l ) calibrations 
involving the N2 (bottom left panel in Fig,ll2[l and the O3N2 
(bottom right panel in Fig. I12p parameters directly reflect 
the ionization structure (see Fig. [7J, having lower oxygen 
abundances [12+log(O/H)~8.10 and 8.05, respectively] in 
the high-ionization region and higher oxygen abundances 
[12+log(O/H)~8.30] in the low-ionization areas. 

The correlation between the ionization degree and 
the results provided by the empirical calibrations is 
shown in Fig. 1131 that plots the reddening-corrected 
[O ill] A5007 / [O 11] A3727 ratio versus the oxygen abun- 
dance derived from empirical calibrations for the spaxels 
that possess an useful measurement of the [O 11] A3727 
emission line. From this figure is evident that empirical 
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Figure 12. Maps of the oxygen abundance deri ved by empirical c alibrations with Ho contours overlaid into the I C 10 PMAS field of: (tot 
left) P ilyugin calibration [average value between iPilvuginl i2001al lbhl using the R23 and P parameters, (top right) Kobulnicky & Kcwlcy 
||2004| 1 calibratio n using the R23 and th e y parameters, (bottom left) cubic calibration between the N2 parameter and the oxygen abun- 
dance deri ved bvlPettini fc Pagell j2004) , and (bottom right) cubic calibration between the O3N2 parameter and the oxygen abundance 
derived bv lPettini fc Page ll J2004) . Symbols are the same that in Fig. [2] 



Table 3. Oxygen abundances, in the form 12+log(0/H), of the regions analyzed in IC 10 using our PMAS data determined using 
empirical calibrations. 



Method 


Parameters 


[HL90] 111 


API (c) 


AP2 (hole) 


AP3 (d) 


AP4 (WR 24) 


Direct 


T e 


8.26±0.09 


8.22±0.08 


8.22±0.17 


8.23±0.09 


8.20±0.12 


P01 


R23, P 


8.15 


8.12 


8.10 


8.13 


8.11 


PT05 


R23, P 


8.13 


8.10 


8.03 


8.13 


8.07 


KK04 


R23, y 


8.25 


8.17 


8.00 


8.22 


8.10 


PP04b 


N 2 


8.22 


8.13 


8.08 


8.18 


8.15 


PP04c 


O3N2 


8.20 


8.11 


8.07 


8.15 


8.12 



calibrations give lower oxygen abundances in regions with 
higher ionization degree, with a monotonic increasing of the 
oxygen abundance w ith decreasing ionization degree in all 
cases except for the IPilvugin fc Thuanl (|2005l ) calibration 
(which changes this tendency for [O ill] A5007 / [O n] A3727 



<1.0). Figu re [131 also s hows t hat the results provided by the 
I Kobulnicky fc Kewlevl l|2004 ) method have a higher disper- 
sion than those obtained using other empirical calibrations, 
and that they tend to be between 0.1 and 0.2 dex higher 
than the others. 
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Figure 13. Relation between the ionization degree, given 
by the reddening-corrected [O III] A5007 / [O 11] A3727 ra- 
tio (x — axis), and the oxygen abundance (y — axis) derived 
from empirical calibrations for the spaxcls that possess a good 
measurement of the [O 11] A3727 emission line. Differ ent sym- 
bols an d colors indicate diff erent empirical calibration s: IPilvuginl 
d2001allbT) (blue crosses) andlPilvugin fc Thuanl l|2005l ) (black A); 
iKobulnickv fc Kewlevl fcooiT using R23 and y parameters (pur- 
ple xl; and IPettini fe Pagell l|2004 l using the N2 (red circles) 
and the O3N2 (green stars) parameters. The oxygen abundance 
determined in the region [HL90] 111 using the direct method, 
12+log(0/H)=8.26, is indicated with a black continuous line. The 
two black discon-tinuous lines shows the uncertainty, ±0.09, of 
this value. 

However, when we integrate the individual spectra of 
an H 11 region, the result provided by the empirical calibra- 
tions are closer to those derived from the T e method. Ta- 
ble [3] lists all oxygen abundances determined using empir- 
ical calibrations in all the regio ns analyz e d here , including 
all the observed area. As we see, IPilvuginl (|2001al lbT) calibra- 
tion provides a constant value of 12±log(0/H)~8.12, that 
is ~0.14 dex lower th an the oxygen abundance determined 
by the direct method. iPilvugin fe Thuanl (|2005t ) calibration 
provides somewhat lower va lues, 12±log(O/H)~8.10. The 
IKobulnickv fe Kewlevl (|2004l l method gives a relatively good 
value for the bright H 11 regions [12+log(O/H)~8.20], but 
the values obtained for the hole and AP 4 (the two spaxels 
where the blue WR bump is detected) are m uch lower. When 
integr ating the flux of all the H 11 region, the lPettiru" fc Pagell 
I2004T ) calibrations using the N2 and O3N2 indices also pro- 
vide oxygen abundances that agree within the errors with 
those determined using the direct method. 

Beside the ionization structure of the H 11 regions, other 
biases (such as the hardness of the ionizing radiation field, 
the structural parameters of the H 11 regions or the star for- 
mation history of the host galaxy) may also affect to the de- 
terminati on of the chemi cal abundances using empirical cali- 
brations (|StasMskall2O10h ■ Hence, we conclude that the anal- 
ysis of star-forming galaxies (and strong starbursting, Wolf- 
Rayet galaxies in particular) using 2D spectroscopy needs 
observations deep enough to detect (at least in the brightest 
H 11 regions) the faint auroral lines, if a good determination 
of the chemical abundances and physical conditions of the 



ioniz ed gas is wanted to be achieved (e.g.. |james et al.ll2009l . 
2010). The deep observations will also provide a better esti- 
mation of the [O 11] A3727 flux, which is essential to precise 
the total oxygen abundance. In the case this doublet cannot 
be observed, deep observations at A7330 are needed to get a 
good determination of the [O 11] AA7318,7330 emission lines, 
which typically are 20-30 times fainter than [O 11] A3727, but 
this spectral range is very affected by telluric lines. Oxygen 
abundance maps obtained using empirical calibrations may 
show features that are not related with the actual metallicity 
distribution of a star-forming galaxy but with the ionization 
structure within its giant H 11 regions. 



6 CONCLUSIONS 

In this paper, we present results from integral field spec- 
troscopy of the intense star- forming region [HL90] 111 within 
the starburst galaxy IC 10. This study is based on PMAS 
data, which cover a wavelength range of 3500-7150 A, with 
an effective spectral resolution of 3.6 A, and map an area 
16" x 16" with a spatial sampling of 1" x 1", that corre- 
sponds to 3.88 pc 2 at the distance to IC 10. We obtain maps 
of different emission lines and analyzed the extinction, phys- 
ical conditions, nature of the ionization, and chemical abun- 
dances of the ionized g well determined locally the age 
of the last star-formation event. 

By defining several apertures, we study the main inte- 
grated properties of some regions within [HL90] 111. For 
the whole FoV, we derive an extinction in the l/-band 
of Ay = 3.05 ±0.11 mag, a total mass of ionized gas 
of (1.76 ± 0.10) x 10 5 Mq, a star-formation rate den- 
sity of 1.64 x 0.09 Mq yr" 1 kpc~ 2 , an oxygen abundance of 
12±log(O/H)=8.26±0.09 and a nitrogen to oxygen ratio of 
log(N/0) = -1.42 ± 0.12. The age of the most recent star- 
formation event is ~3.3 Myr. 

One of our main aims is to precise the places where 
Wolf-Rayet stars are located via the detection of the blue 
and/or red WR bumps. Only two spaxels show an unam- 
biguous detection of the broad He 11 A4686 emission line. 
This feature seems to be produced by a single WNL star, 
its position coincides with [MAC92] 24B. We report a pos- 
sible N and He enrichment at the same location where the 
WR feature is found. Moreover, the enrichment pattern is 
roughly consistent with that expected for the pollution of 
the ejecta of a single or a very small number of WR stars. 
This result suggests that the fresh material released by 
massive stars is quickly incorporated into the surrounding 
ISM. Furthermore, this chemical pollution is very localized 
(~ 2" ~ 7.8 pc) and it should be difficult to detect in star- 
forming galaxies beyond the Local Volume. 

The reddening coefficients derived from the Hcv/H/3 and 
the H7/H/3 ratios have a better agreement when the absorp- 
tion underlying the H I Balmer lines is considered. We find 
a correlation between the electron density of the ionized gas 
and the reddening coefficient, indicating that denser regions 
have higher extinction. Furthermore, the low reddening coef- 
ficients computed around the spaxels where the WR features 
are detected suggest that the gas surrounding the massive 
stars has been evacuated or partially destroyed. In the spaxel 
just to the north of the WR stars location we also detect a 
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broad Ha component underlying the nebular Ha emission 
that may be related to a remnant of a recent SN explosion. 

We discuss the use of the most-common empirical cali- 
brations to estimate the oxygen abundances of the ionized 
gas in nearby galaxies from 2D spectroscopic data. The 
ionization degree of the gas plays an important role when 
applying these empirical methods, as they tend to give lower 
oxygen abundances with increasing ionization degree. We fi- 
nally stress the importance of the detection of the auroral 
lines to derive the electron temperature of the ionized gas 
and obtain a precise map of the chemical abundances in 
extragalactic H n regions. 
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